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Abstract 

We extend the Schwinger boson large iV treatment of the under-screened Kondo model in a 
way that correctly captures the finite elastic phase shift in the singular Fermi liquid. The new 
feature of the approach, is the introduction of a flavor quantum number with K possible values, 
associated with the Schwinger boson representation. The large N limit is taken maintaining the 
ratio k = K/N fixed. This approach differs from previous approaches, in that we do not explicitly 
enforce a constraint on the spin representation of the Schwinger bosons. Instead, the energetics of 
the Kondo model cause the bosonic degrees of freedom to "self assemble" into a ground-state in 
which the spins of K bosons and N — K conduction electrons are antisymmetrically arranged into 
a Kondo singlet. With this device, the large N limit can be taken, in such a way that a fraction 
K/N of the Abrikosov Suhl resonance is immersed inside the Fermi sea. We show how this method 
can be used to model the full energy dependence of the singular Abrikosov Suhl resonance in the 
under-screened Kondo model and the field-dependent magnetization. 
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I. INTRODUCTION 



The work in this paper is motivated by the physics of heavy electron materials. These 
materials are the focus of renewed attention, in part because of the opportunity they present 
to understand the physics of matter near a quantum critical point |jl 0, 0, U, ^|. One of 
the unexplained properties of heavy electron quantum criticality, is that the characteristic 
temperature scale of heavy electron Fermi liquid is driven to zero at the quantum critical 
point □ □□□□ . When either the paramagnet or antiferromagnetic heavy electron phase 
is warmed above this temperature scale, it enters a "non-Fermi liquid" phase. The standard 
"Moriya- Hertz" theory [ll. 12 1 of quantum magnetism is unable to explain the divergence of 



the heavy electron mass in these three dimensional materials. Many other aspects, such as 
the appearance of E /T and H/T scaling in physical properties 



I V o 

14] 



131 . Il4j , the development of 



a quasi-linear resistivity and the tentative observation of a jump in the Hall constant at the 



quantum critical point 



lq suggest that we have not yet found the correct mean-field theory 



for the development of magnetism in these systems. 

These considerations motivate a renewed effort to find the correct mean-field theory that 
spans the quantum critical point between the antiferromagnetically ordered Kondo lattice 



re heav\ 



16, 



12,Ui 



and fully screened Kondo lattice paramagnet. Existing mean-field treatments o 
electron paramagnet describe the spin degree of freedom as a fermionic bilinear 
and while these methods provide an adequate description of the formation of the heavy 
electron bands at low temperatures, they are ill-suited for a description of the antiferromag- 
netically ordered state. 

This suggests that further progress may require a bosonic mean-field description of both 
the Kondo impurity and lattice model. Bosonic spin representations have the advantage that 
they are naturally suited to the description of antiferromagnetism in the Kondo lattice I9J • 
In these approaches, the spin rotation group is generalized from SU(2) to SU(N), providing 
1/N as a small expansion parameter. The hard part of the problem is to capture the 
screening physics of the Kondo effect using the bosonic spin description. A first step in this 
direction was made by Parcollet and Georges . 20], who argued that in order to produce a 
Kondo singlet in an SU(iV) approach, one needs to introduce a multi-channel Kondo model, 
in which the number of screening channels F grows with N. In their approach, it became 
possible to describe the fully screened Kondo singlet by choosing the number of screening 
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channels equal to the number of bosons in the spin representation, F = 2S. One of the 
difficulties encountered in this work, is that the localized moment occupies only 1JN th of 
the singlet, giving rise to a vanishingly small elastic scattering phase shift of ir/N\2u . 

In a more recent return to the Schwinger boson description of the Kondo model 
was shown that a controlled large N treatment of the under-screened Kondo model (UKM) 
can actually be obtained using a single-channel Kondo model. This method captures the 
partial screening from spin S to spin S — 1/2, revealing that UKM is a singular Fermi liquid, 
where the slow logarithmic decoupling of the partially screened moment |24 . 125 ] generates a 
a singular logarithmic dependence of the scattering phase shift [2^ . This is manifested by a 
singular divergence in the specific heat coefficient CyJ T ~ 1 / (T In 4 (T^/T)) and differential 
magnetic susceptibility xi B ) ~ 1/ (5 1^(7^/5)) 0,13 ■ 

One of the interesting aspects of the under-screened Kondo model, is that it displays a 
field tuned Fermi temperature which rises linearly with field. This is a feature found to be 
present at a heavy electron quantum critical point 14]. However, this approach still leads to 
a tt/N phase shift which vanishes in the large iV limit. 

In this paper we continue this earlier work, showing how the phase shift problem is solved 
by introducing K replicas of the Schwinger boson spin, writing 

kN 

where the number of replicas K = kN is scaled with N. This technique preserves a finite 
fraction k of the impurity spin inside the Kondo singlet, and the Abrikosov Suhl resonance 
which develops is now immersed beneath the Fermi sea with an elastic phase shift ttK/N. 
Our results clearly show the formation of a singular Abrikosov Suhl resonance with a finite 
phase shift, but at present they do not extend to the fully screened Kondo model. We shall 
discuss at the end of this article how a future fusion of our method with the multi-channel 
approach may indeed provide a viable Fermi liquid description with a finite phase shift. 

II. THE MODEL 

Our starting point is the Kondo model, which is written 

H = J~] e k c\ a c k(T + JS ■ ipia af3 ipf3. (1) 

ka 



where 5 denotes a spin 5 > |, c\ a creates a conduction electron with wave vector k, spin 
component a, ip a = J2k c ka creates a conduction electron at the impurity site. Our next step 
is to reformulate the UKM as an SU(N) invariant Coqblin Schrieffer model, which enables 
us to carry out a large N expansion of the physics. We write 

H = Yl 6 k C ka C ka + J^Y1 &tf^a ( 2 ) 
ka a/3 

where the spin indices run over N independent values a, [3 £ (1,N). The new feature in 
our treatment, is the introduction of spin replicas. The spin operator is written as a sum of 
K = kN replicas as follows: 

kN 

Saf3 = b lv b fr 

We shall consider the case where there are 25 bosons of each flavor, 

JV 
a=l 

where 5 = 2S/N is kept finite as Af — > oo. The hope behind this approach, is that by 
retaining the bosonic character of the spin, we should later be able to adapt this method to 
describe magnetic behavior in a Kondo lattice. A mult i- channel, rather than a multi-flavor 
formulation of the above model, has previously been treated within an integral equation 
formalism j^j] and a single-flavor version was considered in jsij ] . 

To understand the reasoning behind the introduction of a flavor index, it is helpful to 
consider the strong-coupling limit of this model, where the dispersion of the conduction 
electrons is ignored. When the number of flavors K = 1, only one boson can bind with the 
conduction electrons to form a singlet, and the remaining 25—1 bosons form a decoupled 
local moment, as shown in Fig. ^ (a)- When K exceeds unity, it becomes possible for K 
bosons, each of different flavor, to antisymmetrize and form a singlet with N — K conduction 
electrons, leaving behind K decoupled spins, each with 5* = 5 — 1/2, as shown inU(b). 

The corresponding singlet ground-state is given by 

r N 

i*>= d Vl ...d VK ii(bi^+ci)\s* K ) (3) 

J CT=1 

where we have introduced K Grassman numbers so that the composite operator /J = 
Yli^i^afJIii has the exchange symmetry of a fermion. The ket \S* K ) denotes any state 



N-l 



x 



(b) 



2S-1 




N-K 



V I I 

FIG. 1: Young tableaux which illustrate the strong coupling ground state of (a) the single flavor 
model, with K = 1 and (b) the multi- flavor model with K > 1. In (a), only one boson binds to 
the conduction electrons to form a singlet. In (b), K bosons of different flavor can antisymmetrize 
with each other, to form a singlet with N — K conduction electrons. Weak coupling models flow 
to this strong coupling limit, giving rise to an elastic phase shift 5 = ttK/N. 

formed from 2S* Schwinger bosons of each flavor (e.g. \S* K ) = Ua=i(b\ a ) 2S * |0> )• Carrying 
out the Grassman integral, this becomes 

I*) = eai...a JV &.i-«c&] nw . . -bl KtlK cl K+i . ..cl N \S* K ). 

We shall consider a large N limit in which the ratio k = K/N of flavor to spin degeneracy 
remains fixed, so that for instance, by considering k = 1/3, we make it possible for the 
Schwinger bosons to form a one-third filled singlet ground-state, even in the large N limit. 
The condition that k is fixed as N — > oo guarantees that the impurity Free energy grows as 
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0(N), which is the condition for a controlled large N expansion. 



III. INTEGRAL EQUATIONS FOR THE LARGE N LIMIT. 



Our first step is to cast the partition function as a path integral and factorize the inter- 
action 

JL r i _ . 1 1 _ i 

(4) 



The Lagrangian for this model is then written 



5 = [ dr J2 4 + e *) + £ & ^ + A ) ^ " 2 SK\ + I drH I (5) 

J ° k,a a» Jo 



where the A field imposes the constraint A/& = ^2 afl 2SK. The method we 



now 



follow is closely analogous to that of Parcollet and Georges [20|. First, we integrate out the 
bosons, writing S = S&(A) + Sk, where 



S b {\) = NK \og[l - e~ px ] - 2SKX/T. 



is the part of the action describing the free boson and 



Sk 
Si 



dr 



1 -r 



£ C L ( d r + e*) c k a ~ -jfy^n 



+ Sj 



~m E / drdT^l{r)^{r')D {r - r')0,(r')0 M (r) 



(6) 



is the "Kondo" contribution to the action, where D (iu n ) = iv _ A . 
If we carry out a Hubbard Stratonovich decoupling of Sj, to obtain 
1 rf> 



Si 



PJo 



drdr' 



^V t (T)S(r-r')V(r , ) + E^( r ) II( 



r " r')Mr') 



+ AOI(t-t') 



(7) 



(8) 



D (r - T 

As N becomes large, each term in Sk grows extensively as ~ N, so that in the large N 
limit, the saddle point of this action is expected to saturate the path-integral, giving an 
essentially exact solution to the problem in the large N limit. The Hubbard Stratonovich 
transformation has in essence, replaced 



T — r 



(9) 



These replacements become identities in the large N limit. To see this directly, we integrate 
out the Fermions to obtain the effective Free energy: 



F eff = TS b (X) — A"TTr In 
N 



7 + n 



— A^TTr In [-g~ l + E] 



+ 



P 2 



drdr'nir' - r) 



D a (T-T> 



. T — T 



(10) 



This is the starting point for producing our mean-field equations. If we differentiate w.r.t. 
n(r' — t) we obtain 



KJ(t-t') 

A, 



1 / K 



or 



where 



E(r - t') 



t-t' + 



iV 



A)(r - r' 



-E(t-t') =0 



A" 



iV/3 



^(r-O^r-T') 



i - n M 



Similarly, differentiating w.r.t. E(r' — r) we obtain 

S(r-r') 

r 

1/1 ' (t-t')+II(t,t') 



W 1 " S 



0. 



or 



n (^-^ , ) = -^o(r / -r)^(r-r'), 







where 



1 



£o (w) - E(o;) 



If we convert the integral equations to Matsubara summations, we obtain 

Y,{iu n ) = -kTy^- — -—J(iuj n 
L — ' %v r — X 

r 

U(iuj n ) = ~T^2~ — —g(iu n + iv r ). 



(11) 



(12) 



(13) 



(14) 



(15) 



These equations governing the large A" limit can simply understood diagrammatically as 
"NCA" or non-crossing diagrams, as shown in Fig. |21 
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FIG. 2: Diagramatic representation of large equations, (a) Propagators for the conduction and 
(f) fermion; (b) vertex between particles (c) self-consistent "NC A" diagrams for the self-energies. 

Carrying out the Matsubara summations, we obtain 

E(z) = kn b J( z -X) + k [ ^ 1 ~ f{ - V \ ImJ(y - i6) (16) 

J n z y A 

U(z) = n b g{z + A) + / - nV \j mg{y - iS) (17) 

J 7r z - y + A 

where n& = l/(e /3A — 1) = n(X) and X/A/ = k. These integral equations can be solved simply 
by numerical iteration. 

IV. FREQUENCY DEPENDENT T-MATRIX 

The above integral equations ()16ll7j) were solved by an iterative numerical procedure. 
One of the key quantities of interest, is the conduction electron phase shift, given by 

S c = Imln [1 — 27rp£(0 — iS)] , 
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Under the assumption that K = kN bosons are bound into the singlet, Friedel's sum rule 
determines that the conduction electron phase shift will be 5 C = kir. Our numerical results 
(Fig. EJ) confirm this result. 
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FIG. 3: Showing the numerically computed dependence of the conduction electron phase shift on 
k. These phase shift were actually computed from at zero temperature from the zero-field limit of 
the finite temperature conduction electron Green's functions. (See section VI) 



To get an approximate understanding of the numerical results, it sufficient to carry out 



the first iteration. If we assume g 



(0), 



CO 



iS) = inp, where p is the conduction density of 



states, then the first order approximation for IT(z/) is 
H^(u — i5) = iTcprib + pin 



pj{y - iS) 



P 



v — iS 

./,, n (^-^) \ u (Tk\ + i n [ nb + 0(- v )] 



(18) 



corresponding to an effective Kondo interaction J{y) which changes sign, from positive and 
antiferromagnetic (Re[J] > 0) at high energies to negative and ferromagnetic (Re[J] < 0) 
below v ~ Tk- Since the interactions become weak at both low and high energy, this 
expression captures the essential character of the full solution in these limits. To fine- 
tune the solution, we do however need to take account of the renormalization of the local 
conduction electron density of states at low energy. The renormalized conduction electron 
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propagator at the Fermi energy is given by 

go rip 



9(0 - i5) 



l-g o E(0-iS) l-inpE(p-iS) 

Now we can determine S(0) by noting that at the Fermi energy, the scattering is elastic, so 
that S(0) is real. From the conduction electron phase shift, 



it then follows that 



5 C = irk = Imln [1 — 27rpX(0 — i8)] , 
tan(7r/c) 



S(0) = - 



Tip 



from which we deduce that 



1 - i7rpS(0 - i8) = e i7rk sec(7ik), 

g(0-i6) = i7r pe~ ink cos nk (19) 

so the renormalized density of states is given by p* = ^lm[g(0 — i5)} = pcos 2 (irk). This 
depression in the local density of states will reduce the coefficient of ln(z/) in the frequency 
dependence of the inverse coupling constant. If we approximating the renormalized density 
of states by p* within an energy Tk of the Fermi energy and p otherwise, we may write an 
improved approximation for II (z/), as 

T K \ D 



n (2) (z/-i<5) = p* ( inn b + In ] +pln 
\ v — id J 



pJiy - is) = l p - = - 7 , sec2(7rfc) , (20) 

_ _ n(v - id) ln ^ + m[nb + 9 (_ v y 

The logarithmic scaling and the dependence —j ~ cos 2 (7rfc) ln(T^-/|z/|) are indeed confirmed 
in our numerical results (Fig. 0}. One of the interesting points about this result, is that 
the propagator is logarithmically dependent on energy and does not develop a power-law 
dependence on energy that is the hallmark of the overscreened Kondo model[22jL and also 
appears to be present in a recent fermionic parallel to the current approach [23|. The x 
fermion basically represents the singlet combination of conduction electron and Schwinger 
boson Xn ~ The absence of a powerlaw in the x propagator is evidence that 

the addition, or removal of a boson from the system does not change the scattering phase 
shift. This is presumeably because the addition or removal of bosons from the system merely 
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FIG. 4: Showing logarithmic dependence of Re[l/J(0)} on magnetic field at zero temperature, 
calculated numerically for the case k = 0.2. The magnetic field provides the cut-off to the loga- 
rithmic scaling. For details of how the magnetic field was introduced, see section VI. Dashed line 
is the the curve y = cos 2 (irk)ln(^-) . 

changes the size of the undescreened moment, without altering the number of bosons that 
are bound into the singlet. 

At low temperatures J scales to zero with logarithmic slowness. If we look at the the 
conduction electron self-energy as given by (JIEJ), we see that the leading singular frequency 
dependence is given by the kn^JiiS) term, so that at low frequencies, 

k(nf, + 9{v)) sec 2 (nk) t&n(nk) tan(5(a;)) 



EM 



ln( St 



i7r[n b + 0(u)} 



np 



np 



(21) 



where 



5 c (uj) ~ nk 



1 + 



+ o 



\n\T K /\uj\ 



(22) 



H t k/H). 

This singular energy dependence of the conduction electron phase shift is a consequence of 
coupling between the degenerate manifold of the partially screened moment and the con- 
duction sea. Although the scattering of conduction electrons is elastic at low energies, the 
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singular frequency dependence displayed here sets this system apart from a conventional 
Landau Fermi liquid (where the phase shift depends linearly on energy). Indeed, if we use 
E(u;) to define a frequency dependent wave-function renormalization constant, we find that 
this quantity diverges as Z = j^tt^ ~ [Jk^yAjj > so that there are no well-defined quasi- 
particles associated with the Kondo scattering. This is the meaning of the term "singular" 
Fermi liquid. 

The singular energy dependence of the scattering is also reflected in the shape of the 
Abrikosov Suhl resonance, given by the imaginary part of the electron t-matrix 

Im[t(u)] = Im 

These basic features are each borne out in the detailed numerical solution of the integral 
mean-field equations of the large N limit. Fig. El shows the results of this numerical 
calculation, showing the spectral function of the t-matrix and the dependence of the phase 
shift on the number of flavors. 

Im t (<y) 



71 p Im t (a) ) 

sin 2 (7E k) o*s 

0.4 
0.2 

-7.5 -5 -2.5 2.5 5 7.5 10 

a>/T K 

FIG. 5: Showing the frequency dependence of the t-matrix, normalized with respect to zero 
temperature value at Fermi energy. In this numerical calculation, k = K/N = 0.45, n& = 2S/N = 
0.2. Numerical labels give temperature T/Tk- 

12 



1 - g (u)E(uj) 



Tip 



sin 2 (7r&;) + 



irk[rib + {(*))] sin(27r/c) 
ln(7WM) . 



(23) 




This singular energy dependence of the t-matrix also manifests itself in the temperature 
dependent resistivity, given by 

P(T) - mirp f duj ( -|£ ] Imt{u) (24) 



Pu J \ OUJ 

where pu = (^~^J is th e unitary resistivity as shown in Fig. El 




V. PHASE SHIFT AND THE SCREENING OF THE MOMENT 

To reveal the screening of the local moment, we need to compute the first correction to 
the constraint equation. To satisfy the constraint we must differentiate the Free energy |T0|) 
w.r.t A. This yields 

2SK = KNn b + 5M b 

5Af b = —T 2 K c _^ 2 g{ivn + iu r )J{iu r )- (25) 

\%V n A) 

The first term in this sum determines the number of bosons that condense into the unscreened 
magnetic moment. The second term 8Mb is the total number of bosons bound into the Kondo 
singlet in the ground-state. 

13 



At zero temperature we can replace the discrete Matsubara sums by a continuous integral 
as follows, 

^ I ■ 2m ^ .L^2m 

so that the number of bound-bosons becomes 



where we have replaced ^ u \y2 —> (^hyj^j- Next we can integrate the internal integral 
by parts 

AA4 = -K f -gL f *L * {g{y + u)) J-J( U ), (27) 
J (2m) J (2m) du v — A 

and then replace 4- — > 4- to obtain 

1 du dui 



n(o 



W > = K J^l{-J&)« v + a) ^)™- <28) 

We can now identify the term inside the central brackets as the self-energy IT(u;), which 
enables us to compactly rewrite the integral as 

= -- [fy(0) - fy(-oo)] = -—S* (29) 

7T 7T 



where 



so that 



5^(00) = Imln 



- - U(u - i5) 



2KS* 

K2S = / NKr? b <L . 

7T 

Now the argument inside the logarithm that determines 5<f, is the inverse of the running 
coupling constant, 

-J— = \-U{w-i$). 
J{oj) J V ; 

In the underscreened Kondo model, the residual coupling between the partially screened local 
moment and the conduction electrons becomes ferromagnetic at low temperatures, scaling 
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logarithmically to zero. This means that j^_ iS ^ ~ — P 1 hi(T^/u;) — iS at low energies, which 
implies that 5$ = —n. This in turn, implies that the number of bound-bosons is 5Mb — K. 
As a consequence the number of unpaired bosons is reduced by K in the ground-state, and 
if we define S* = y^, the effective spin of each flavor will be 

We can also relate the number of bound bosons to the conduction electron phase shift 5 C , in 
a parallel fashion, as follows: 

-E(w) 

, A , 

du d ( f dis „. , 1 



„ /■ du SEfw) , . 



= - [<J C (0) - <Jc(-oo)] = -5 e . (30) 

7T 7T 

In deriving the third line of this expression, we have taken the large band- width limit, 
enabling all frequency dependence of the conduction green function g a = —iirpsign(uj n ) to 
be be ignored. The conduction electron phase shift is given by 

8 c {oj) = Imln [1 — g (uj — z5)S(u; — i5)] . (31) 

By comparing the two expressions (J2o1 ) and (j3*Ul ) for 5Mb, we are able to confirm that the 
conduction electron phase shift is given by 

,5M b \ 

7T — — = 1XK. 



N 

Notice that along the way we have proven that 

N (I) + K ( » ) = (32) 

Although we have proven this strictly in the large N limit, this result is a Ward Identity that 
is expected to hold for all N, a result which relies solely on Fermion number conservation 
(see Appendix A). This result is in effect, a statement of famous the "Anderson- Clogston 
compensation theorem" - that the total number of bound fermions bound by the Kondo 

r 

effect in the infinite band-width limit, is zero| 
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VI. ZERO TEMPERATURE MAGNETIZATION 



In order to examine the effect of a magnet field in the large N limit we need to be careful 
about our definition of the magnetization. We shall suppose that the magnetic field couples 
preferentially to the "up" spin (s) of each flavor, i.e. that the magnetization takes the form 

M = gN{n b ^ - n b ^\ B =o) 

There are various ways in which we can now definite n^. One way to do this, is to define 
the "up" states as the first K spin components, i.e 

K 



M« = 



N 



2(N -K) 



2SK 2 



N 



Here, value of g n is chosen so that at full polarization, M = SK. With this definition, 
the first K spin channels of the conduction electron are "up" electrons and the remaining 
N — K spin channels are "down" . When we add to the corresponding expression 

for the conduction electrons, we form a conserved quantity, so that we shall call the 
"conserved magnetization". By imposing the constraint, XIo-^^^^m = ^SK, the conserved 
magnetization can be rewritten as 



M« = - 
2 



2SK 



N 



N-K 



where n^i = 2SK — = X]«=i 2~2 a >K n ba^- Now unfortunately, if we couple the magnetic 
field up to then we do not completely remove the spin degeneracy of the ground-state. 
For this purpose, we need a more restrictive definition of n^, and we shall choose "up" to 
mean the spin component where fj, — a, i.e 

K N K 

H=l cr=l /I=l 

With this definition, by imposing the constraint, J2 all n ba^ = 2SK, we obtain 



M 



We shall use this definition as our method for coupling the magnetic field to the spin, so 
that M = —5F/5H is the true "thermodynamic" magnetization. In the case where K — 1, 
the thermodynamic and conserved magnetizations correspond exactly. However, the need 
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to break the flavor or replica symmetry at finite K unfortunately forces us to delineate 
between these two forms. However, we shall see shortly that (M^) and (M) have almost 
identical expectation values, and that the conserved magnetization has a far more convenient 
expression in terms of the scattering phase shift. The coupling to a magnetic field then gives 
the Hamiltonian 

H B = -BM = -BKS + ( 33 ) 

Now in the large N limit at low temperatures in a field, the "up" bosons condense. This 
allows us to carry out the Kondo version of spin-wave theory. Formally, we condense the 
"up" bosons and integrate over their phase fluctuations to exactly impose the constraint. 
The resulting Holstein Primakoff transformation is obtained by replacing 



bm _ < (' = «)• (34) 

b aa (a ^ a) 

where rib a = ^2 a ^ a ^aa^ a - 111 this process, the a = a or "up" components of the boson 
fields have been eliminated. In a field, the amount of fluctuations n^a ~ 0(1), so that as in 
spin-wave theory, to leading order, we can drop the n& a inside the square-root. With this 
understanding, inside the interaction, we must now make the replacement 

P 

so that the Hamiltonian in a field now becomes 
H = ^ e k c\ a c ka + Hj + Hb 

ka 



J 



(35) 



The effective action in a field now becomes 

1 _ 1 

Si + S n 



I 



P 

dr 



£ 4 ( d r + e fe) C ^ ~ " SK(2X + B) 



k.a 

rP 



Si = ~^££ f ' drdr^i(T)Mr')D a (T - r')i(r')i(r) 

a=l u+a ^° 
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I2S K f 13 

S mix = dT ^(t)10(t) q + 4> a {r)iP a {T)] 

a=l J° 



(36) 



where D a {iv n ) 



The interaction term Si can be split up into a term with no 



iu n -(B/2) ■ 

restrictions on the spin flavor summations, plus a term 1Z that can be neglected in the large 
N limit: 

K N r p 



a=l a=l ^° 

U = I drdT^i(T)Ur')DAr-r')Mr')Mr). 

a=l ^° 

If we carry out a Hubbard Stratonovich decoupling on the first term, we now have: 



(37) 



S= I dTdr'C(T,r') - (3SKB 
Jo 



K 



C = [-go 1 + E] ^ + E W' 



ok 



a=l 



-# 1 + £ v 7 ^ 

- j- 1 + n 



'n 6 



+ Nn(T-T')———X(T'-T) 



(38) 

(39) 



^o(r-r') 

where fi h = 2S/N Taking the saddle point values for n and E, and then integrating out the 
fermions, the mean-field free energy in a field is then 



NT 



= -(1 - A;)Trln [-g' 1 + E] - A;Trln 



-9o 1 + E 

-j _1 + n 

+ N f dTdT'n(r',T)-—^ -E(r,r')-^ 5 - 

7 jJo{t - r) 



(40) 



When we impose the saddle-point condition on E and n, decomposition, we must be care- 
ful to delineate between the propagators of "up" and "down" electrons. The mean-field 
equations are then 



E(r-r') = -kD(r - t')J(t - r') 

n(r-r') = -7J(r'-r)[% T (r-r') + (l-% i ]. 



(41) 



The mixing terms between the "up" electrons and the "phi" fields mean that we must now 
modify the propagators as follows: 

0i M 



9o V) - E|H 



18 



Here 



(43) 



and 



J(c) 



J 1 - n(w) - n 6 ^ (w) 

where h b = 2S/N. 

If we convert the integral equations to Matsubara summations, we obtain 
£(^„) = -* £ iVr -\ B /2) J{iUn - Wr) 

Carrying out the Matsubara summations, we then obtain 

i - /(y) 



(44) 



(45) 



E{z) = k J 
U(z) = 



7T z - y - (5/2) 

dy f(y) 



ImJ(y — id) 



Im [kg^y - iS) + (1 - - i<J)] 



(46) 
(47) 



vr z-y + (B/2Y 

From the solutions of these equations, we can compute the field dependent propagators 
and free energy. We can define the following scattering phase shifts associated with the 
conduction electrons and fermion: 



84(B) = Im 



5 ca (B) = Im 



ln(l- J[U(z)+n m (z)\) 



z=0-i& 



ln(l -goE ff (z)) 



("=T, I) 



(48) 



z=0-iS 



These phase shifts at finite field are actually related by Ward identities (see appendix B), 
and obey the following relationships at all fields: 



K5 4> + N5 [ = 
K5j + (N — K)8 i = Ktt 



(49) 
(50) 
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The differentiation of the Free energy w.r.t. field is identical to the differentiation w.r.t. 
A carried out in the previous section. The magnetization is given by 

M = SK - 

2 

5M h = —T 2 K (■ _„ /0 v 9(wn + iUr)J(iUr), (51) 

where g = kg^ + (1 — k)gi- If, instead of using the thermodynamic magnetization, we use 
the conserved magnetization then we find that 

M (1) = SK 



SAf b 



2 

5Af b = -T 2 K n _ B/2 )2 9l( iU n + iu r )J(iUr), (52) 

The difference between the two expressions results from the slight energies. The advantage 
of the second quantity, is that by being related to a conserved quantity, it is related to 
the scattering phase shifts. Using exactly the same methods we relate &A/& to the 4> and 
conduction electron phase shifts, as before, to obtain 



KS^B) 

KS- 



M (1) = KS + 

2 7T 

NS^B) 



= K {3 -\)^m^m ( 5 3) 

where we have used the phase-shift identities ()49)1 given above. 

We have solved equations f)42l E3*l at zero temperature, finite field B by numerical 
iteration, using a Fast Fourier transfer routine to carry out the convolutions. Fig. [7| shows 
the field dependence of the coupling constant g(B) = Re[J'(uj)\ u , = o, showing the change in 
sign of the coupling constant as the system goes from weak coupling at high fields to strong 
coupling at low fields. Fig. El shows the field dependence of the "up" and "down" phase 
shifts. The results obtained by our large iV method are strikingly similar to results recently 
obtained by numerical renormalization group and the Bethe ansatzj^. Fig. shows the 
field dependent reduction in the magnetization AM = —K8 ( j ) (B)/7v. 



VII. DISCUSSION 



Our method gives a controlled treatment of the underscreened Kondo model, and using 
the method of replicas, we have been able to develop a treatment of the Kondo effect which 
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FIG. 8: Field dependence of the "up" and "down" phase shifts for k = 0.45. 

correctly reproduces the finite phase shift 5 = nK/N produced by the scattering resonance. 
What are the prospects for going further? 

The first question of interest, is whether this approach can be used to describe the fully 
quenched Fermi liquid? An adhoc way to attempt to access the Fermi liquid, is to seek 
solutions where 

M b = K = 2KSn b {\) + N^- 

7T 

The presence of the additional term in the constraint then raises the value of A, so that in 
the ground-state A ^ 0, and hence rib 0, producing a fully quenched ground-state. The 



21 




Log w (B/r K ) 

FIG. 9: Field dependence of the magnetization reduction AM/K = d^/ir for k = 0.45 

obvious drawback with this approach, is that by pushing the problem to this limit, we are 
exploring a limit where the occupancy of each boson spin/flavor state is of order 0(1/N), 
which lies outside the strict region of validity for the large N expansion. 

A more encouraging approach to the problem may lie in trying to unify the boson replica 
approach used here, with the multi-channel approach developed by Parcollet and Georges 
(PG). The PG approach introduces F flavors of conduction electron, considering an interac- 
tion of the form 

j F 

R i = (54) 

7=1 a(3 



(a) 



(b) 



-fN 



t 

kN 
I 



N 



2S= fN 



-fN- 



FIG. 10: (a) Rectangular Young Tableau which will self-assemble in strong-coupling ground-state 
(b) of color-flavor-spin model in the case where F = 2S. 
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This approach leads to a scattering phase shift 5 = tt/N which vanishes in the large N 
limit, however, it has the virtues that for F = 2S it does describe a fully quenched Fermi 
liquid, and moreover, the method can, in the lattice, be neatly combined with the Arovas- 
Auerbach jiol l tr eatment of antiferromagnetism, using a SP (2N) description of the RKKY 
30, 



interactions 



Is it possible to extend this model, introducing both conduction flavor and boson color 
at the same time? This suggests the following color-flavor-spin (CFS) model 

H ' = a7 E (^)W e M. V G l^K]) (55) 

where we choose 2S = F = f'N, and K = kN in the large N limit. The strong-coupling 
solution to this model involves the formation of a singlet between a rectangular representa- 
tion of the Bosons, and the F channels of conduction electron, as shown in Fig. El in the 
ground-state, we expect the bosons to self-assemble into this strong-coupling representation. 

To develop a large N expansion of the CFS model, we might consider writing down the 
connected skeleton graphs that enter into the Luttinger- Ward functional Y [G, <J, D\. 
The two leading diagrams in Y take the form shown in Fig. ^1 There is one loop for each 



Y(G,J,D)= 





OCN 2 ) O(l) 



FIG. 11: Two leading order skeleton diagrams for the Luttinger Ward functional of the CFS 
model. 

quantum number in each of the above diagrams, but because the second diagram involves 
more vertices, it is smaller by a factor of 0(1/N 2 ). Y can be used to generate the self-energies 
for a conserving approximation, via the relationships: 

5Y 



NFE(iu n ), 
oG{iu n ) 

5Y 

KFU(iuj n ), 



SJ(iuj n ) 
SY 

—— = NKZ b (iis n ), (56) 
6D{w n ) 
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where the pre-factors arise because of the sum over spin, flavor and channel that arises 
inside Y. The first term in this series can be used to generate the diagrams that interpolate 
between the PG and the replica approach. This term in the skeleton expansion of the Free 
energy is of order 0(N 2 ) in an approach that involves color, flavor and spin. Since next 
skeleton diagram is of order 0(1), so that at first sight, we may neglect all higher order 
diagrams. 

However, it turns out that the development of a strict large N expansion encounters a 
technical difficulty with the profusion of higher order planar diagrams. Similar difficulties 
have been encountered in the large N treatment of quantum chromodynamics (QCD). To 
see this, it is useful to relabel the conduction, boson and <fi propagators as two parallel lines, 
carrying the respective quantum numbers of spin, color and flavor, as follows 



a, a o 

— > = > 

a 




a 



With this notation, we see that introduction of the additional flavor quantum number, 
means that the interaction vertex between bosons is enhanced by the sum over virtual flavor 
fluctuations, from from a 0(1/N 2 ) term in the current theories, to a 0(1/ N) vertex as shown 
below. 




oa/N) 

FIG. 12: Magnetic vertex. The additional flavor index enhances this vertex by a factor of N. 
If we now look at the boson self-energy diagram generated by this vertex, 
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we see it is of order 0(1). Furthermore, when we close the external legs on this diagram to 
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create a skeleton graph, 




0{N 6 /N 4 ) = 0{N 2 ) 



we see that this graph (which contains six internal quantum number loops, ) is of order 
0(N 2 )- the same order as the leading graph. Unfortunately, as in QCD, this is just the 
beginning of an entire profusion of higher order "planar diagrams" which are all of the same 
order in the large N expansion. 

Despite this difficulty, it may be that the leading order 0(N 2 ) diagram shown in Fig. 9 
is already sufficient to generate a good conserving mean-field theory for the fully quenched 
Kondo model. We do not yet know how the neglect of the higher order planar diagrams 
affects the results, but it is tempting to speculate that these diagrams are irrelevant in the 
Fermi liquid, or magnetic ground-state, since the flavor and color quantum numbers become 
massive when the Kondo singlet forms, or when there is an applied magnetic field. In this 
case, the higher order diagrams may only renormalize the leading order term in the skeleton 
graph expansion of the Free energy. These considerations lead us to suggest that the leading 
order skeleton free energy diagram for the color-flavor-spin model may provide the key to a 
successful mean- field theory that spans the magnetic quantum critical point. 

One of the interesting final questions that deserves discussion, concerns the physical 
significance of the additional quantum numbers that we have introduced. The color-flavor- 
spin large N approach is basically approximating the single box of the Young Tableau 
representing a spin 1/2 by a rectangular Young tableau. What is the meaning of the intensive 
variables / = F/N and k = K/N that appear in the PG and the current approach? One 
fascinating possibility, is that quantum numbers associated with these degrees of freedom 
describe the internal quantum numbers of the composite quasiparticle in the heavy electron 
state. In the Fermi liquid, we expect these quantum numbers to be inert, but at the quantum 
critical point, these quantum numbers become unconfined. These issues will be followed in 
forthcoming work. 

This research is supported by the National Science Foundation grant NSF DMR 0312495. 
We should particularly like to thank Anirvan Sengupta for discussions concerning the matrix 
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aspects of the CFS model, and the profusion of planar diagrams. Shortly after submitting 
this paper to the archive, we became aware of a closely related fermionic "replica" approach 
by S. Florens, which we have referenced in this revised draft. (^J Discussions with N. Andrei, 
G. Kotliar, P. Mehta, G. Zarand and L. Borda are gratefully acknowledged. 



VIII. APPENDIX A 



In this appendix, we derive the Ward Identity 

<§) + *©=° 

and relate it to the Anderson- Clogston compensation theorem. According to this theorem 
the net polarization of the conduction sea by a localized impurity is zero, in the limit 
of a broad band. A localized resonant scattering center induces Friedel oscillations in charge 
density of the medium, and tends to reduce the charge density in the immediate vicinity of 
the impurity. The compensation guarantees that this local depression in charge density is 
compensated by an enhancement at greater distances. 

To understand this compensation effect, we need to examine the conduction electron 
Greens function, which is given by 

G k) k>(iuj n ) = h t[iu n )- 

iu n - e k iu n - e k iu n - e k / 

where t(iu n ) is the t-matrix of the impurity. The second term in this expression induces the 
Friedel oscillations in charge density, which are given by 

Ap(x) =NT Y t(iu n ) e <(k-k')-x e ^ l0 + 

. tri u*> n - e k iu n - e k / 

K,K ,lU n 

For example, the change in density at the impurity is given by 

Ap(0) = -A(ttp) 2 / —f(u)1mt(u - i8) < 
J-oo n 

( where we have carried out the momentum sum, replacing ^ k (^„ — ^k) ~^ —iirpsgn(n)), 
corresponding to a reduction in the electron density. By contrast, the total change in density, 

AQ= f d 3 xAp(x) = NT Y -A{iuj n y^ (57) 

K.,lU) n 
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Assuming that the width of the resonance is much narrower than the bandwidth D, we can 
replace 



D 



T vl^P de T ^ = 7T + °"m 

(^n-e k ) 2 J-D (iuo n -tf D D 2 



so that 



A Q = ^^tK) e - 0+ = -L „ o (58) 



iui n 



where 



2np f — f(u)lmt(u — iS) 

J-oc 7T 



is a measure of the width of the resonance. Thus the total polarization of the conduction 
band is of order Y/D « 1, which is negligible in the limit of infinite bandwidth. 

Now we shall relate AQ to the scattering phase shifts. To do this, we appeal to the 
Luttinger-Ward functional Y[G, J\lor this problem, represented by a sum over all skeleton 
closed loop diagram contributions. The differential of this functional with respect to the 
Greens functions generates the corresponding self-energies: 

- t—t. — r = £<r 
bg a {iuo n ) 

5Y 

- YTT- \ = U ^ luJ n), (59) 

where we have explicitly displayed the spin and flavor indices. Now the conservation of 
charge guarantees that each of these can be decomposed into one or more closed fermion 
lines. In the low temperature limit, the Matsubara sums along these lines may be replaced 
by continuous integrals along the imaginary axis, 



lLO n 



dz 



Now if, at zero temperature, the frequency running along each such loop is incremented by 
a small amount z — > z + i5uo, Y is unchanged, because the the shift in z can be absorbed by 
a simple change of variable. We deduce that at T = 0, in the absence of a field, 



5Y = i5uo 



/duo 
2m 

= —ibuo [ - — : ( 
J 2m V 



SY dg a 5Y dj^ 

5g a (oo) duo dj^uo) duo 

N^iuo)^- + KU(oo)^f- J = (60) 
duo duo J 
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Integrating this result by parts, we obtain 

= 



This valuable result is a consequence of fermion, or charge conservation. 
Now the total change in the charge of the system is given by 

/i 

lU>n 

where the final integral is along the imaginary axis. 

We now subtract the result ()61|) from this expression, to obtain 



(61) 



•o+ 



(62) 



= N / — 



j0+ 



T ^ f du .dll , 
K / —J(u>)—& 
J Ziri duo 



•0+ 



N 



N 



I til 0»P - " 11 ■»!-*-']) ^ + Kf±lH-J-\W 

duj d 



o+ 



27ri da; ' y 2ni du 

Finally, we distort the contour around the negative imaginary axis, to obtain 



(63) 



AQ = / — Im— 

7T du 

X ( ^ ] + K I ^ 

7T / V 7T 



JVln[l - jo(z)E(z)] + ATln[l - J U(z)} 



-id 



(64) 



where the phase shifts are defined as 



6 C 
8a 



Im 
Im 



In (1 — iTcpL{uj — i8)) 
ln(l - J TI(uj-i5)) 



uj=0 



(65) 



oj=0 



Using the compensation theorem to set AQ = in the infinite band width limit, we obtain 
the sum rule 

= NI S A + k(*±) (66) 



IX. APPENDIX B 



In this section, we prove the finite field Ward Identities, 

K8 <p + N5 l = 
K5 } + (N- K)5 l = Kir 



(67) 
(68) 
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The first of these results is a finite field generalization of the result of Appendix A. In a 
field, the "up" conduction electrons and the (f> fermion become hybridized, so that and 
J are no longer independent variables. It is convenient to introduce 

1 



Jo 1 



nM 



in terms of which 



J(u>) 



1 



J -n h gi{u) 
1 



(69) 
(70) 



h b J{u) 

When we write the Luttinger Ward functional, we must be careful to express it as a function 
of independent propagators. One possible choice is g± and J . In this case, we can divide 
the continuous fermion line running through Y into sections which are either g± or J . (In 
this procedure we have effectively integrated out the conduction electrons first, so that the 
fermion propagator contains a contribution from its hybridization with the conduction 
electrons in a field. ) Sections that involve the "up" electron propagator can be broken up 
into gi and J using (|70|). If we then write Y = Y[gi,J~] and vary the frequency running 
along the continuous fermion line, we obtain the finite field version of (JETI1 ). 



SY[g u J] 



-i5oj 



2irt \ duo duj 







(71) 



Following the same steps that were taken in Appendix A, we obtain the first of relations 

(EH), 



where in a field, 



SJB) 



Im 



Im 



K5^ + N8 l = 0. 



ln(l- J[U(z)+n m (z)\) 



=0-i<5 



ln(l-SbSi^)) 



(72) 



=0-i<5 



Now alternatively, we can take the independent propagators to be g^, g^ and J . That is 
to say, we are effectively first integrating out the fermion, so that its propagator J does 
not contain a contribution from hybridization with the "up" conduction electrons, whereas 
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the " up" electron lines now contain a contribution due to hybridization with the fermions. 
When we make a variation of the frequency along fermion lines inside Y, we obtain 



dY[g h g h J] = ibuo 



duo 
2tH 



o, 



where E| = Ej + UbJ, so that 
= - 



duo 

2?7i 
duj d 

2ni du 



Kg,(uo)^l + (N- K) 9i {u) d ^ + KJ{uo) 



-S 



du 



du duo 
K In (1 - g ^{uo)) + (N-K) In (1 - # £|M) + Kin (1 - J U(u)) 



71 







(73) 



where we have defined 

8$ = Imln [1 - J U(z)] \ z= - iS 

Now since U(—i8) = ReU(0) + i5, 8^ is or — ir depending on the sign of J(0). We may 
write 

K8^ + (N — K)8 i = K(n + s) 

where s = |(1 + sgn( t 7(0))) = 7mln(— t 7(^))| z =-i < j. In actual fact, 5f will jump by 7r exactly 
at the point where J changes sign, so by redefining 

1 — inpT.^z) 



5t = Imln 



we obtain a phase shift that evolves smoothly with field, which satisfies the Ward identity 



K5-[ + (N — K)8 l = Ktc. 
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